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Abstract 

A path integral formalism for multispecies anyons is introduced, whereby partition 
functions are expressed in terms of generating functions of winding number probability 
distributions. In a certain approximation, the equation of state for exclusion statistics 
follows. By Monte Carlo simulation, third-order cluster and virial coefficients are found 
numerically. 

1 Introduction 

In two dimensions, a concept of statistics of distinguishable particles (multispecies statistics, 
mutual statistics) arises ^ , by close analogy with the concept of fractional (anyonic) statistics 
12]. Multispecies anyons are particles whose wave function picks up a statistical phase when 
a particle of one species encircles one of another species. Multispecies statistics is interesting 
per se, as a novel basic notion in quantum mechanics, as well as in view of its application to 
the quantum Hall effect, because quasiparticles and quasiholes in quantum Hall systems do 
in fact obey such statistics jHj. 

A particular model of multispecies anyons in a high magnetic field has been solved exactly 
I^IHj. The generic case of free multispecies anyons, however, has not been explored beyond 
the (trivial) 2-body system. This is what makes the subject of the present paper, the focus 
being on the statistical mechanics of the free anyon gas. 

For the system of more than two identical anyons, there appears to be no exact quantum- 
mechanical solution (multiparticle energy levels cannot be expressed in terms of single-particle 
ones like they can be for bosons and fermions), and even the high-temperature cluster and 
virial expansions from the third order onwards are only known approximately. One cannot 
expect to do any better with multispecies anyons. Two distinct approaches have been suc- 
cessfully employed in the treatment of identical anyons: direct numerical computation of the 
spectrum JHIIZI and Monte Carlo calculation of partition functions JHIISI- It is the latter that 
we resort to here. 

After introducing the basic definitions and relations in Sec. 2, we develop, in Sec. 3, 
a representation of partition functions in terms of generating functions, which are Fourier 
transforms of probability distributions of winding numbers. This is a direct generalization 
of the technique developed in Refs. |H1-Cni- Then, in Sec. 4, we introduce an approximation 
jH] in which only two-particle correlations are accounted for. It yields simple polynomial 
expressions for the generating functions, resulting in a virial expansion which depends on 
the statistics parameters at the second order only. This corresponds to the equation of state 
for multispecies exclusion statistics (21,^5 • Section 5 contains the results of Monte Carlo 
simulations of the generating functions and the corresponding mixed virial coefficients. 
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2 Definitions and basic relations 

Multispecies anyons are particles in two dimensions, characterized by a statistics matrix 
||aafe||, such that interchanging two particles of species a supplies the wave function with a 
phase factor of exp{nTaaa) and pulling a particle a around a particle b supplies exp(2i7ra(ib). 
Due to translation invariance, the matrix is symmetric. Due to periodicity, it is enough to 
consider < aaa < 2 and < Oab < 1- 

The equation of state of a multispecies gas is obtained as a straightforward generalization 
of the single-species case. The fugacity expansion of the grand canonical partition function is 

H= E Z^,...^,z^---zf, (2.1) 

Ni...Ns 

where .Z'7Vi...Afs is the canonical partition function of the system containing Na particles of 
species a for each a = 1, . . . , s, and Za = exp(/3/ia) are the fugacities (no chemical equilibrium 
between species assumed). The cluster expansion is (P is the pressure, A the area) 

/3^=^= E bM,...N.z^^ ■ ■ ■ z^' (2.2) 

Ni...Ns 

and the virial expansion is 

f5P= Y. Ak,...k.Pi'---Ps\ (2.3) 

where pa = Zad{f3P)/dza are the partial densities. Relations between partition functions, 
cluster and virial coefficients are easily deduced. The cluster coefficients up to the third order 
are 

Ab2 = -^^1 + ^2 , ^&ii = -^01^10 + Zn ; 

^^3 ~ o^i ~ ^1^2 + Zs , Ab2i = ZqiZ^q — ZioZii — Z01Z20 + Z21 , (2.4) 

Abiu = 2ZqqiZqiqZiqq — ZiqqZqii — ZqiqZiqi — ZqqiZuq + Ziii , 

and the virial coefficients are 

A2 = -B2 , All = -Bii ; 

A3 = 4P2 _ 2B3 , .421 = Pn + 4^11^20 - 2^21 , (2.5) 

^111 = 2i?oii-Bioi + 2i3oii-Biio + 2i?ioi-Biio — 2i?iii , 



where 



Bm...Ns = ,^1 .Z''"''^ .iv. ■ (2-6) 

"lOO...O"010...0 ' ■ '"O.-.O! 

Dimensionless virial coefficients are 

"-ki-ks = ^2{ki+-+L-l) ' ^^-"^^ 

A = \/2'ir(3/m being the thermal wavelength. (Cf. [Hini; the cluster coefficients here contain 
an extra factor of 1/A as compared to the notation used there.) 



3 The generating functions and their connected parts 

We now wish to represent the partition function Zni...Ns as a path integral over closed paths 
in configuration space. Since particles of the same species are indistinguishable, a path may 
induce a permutation of such particles. This breaks up all paths into classes which constitute 
a representation of the "colored permutation group" , a direct product of permutation groups 
Sni fS) . . . <8> Sns- To a conjugation class V of elements of this group there corresponds a set 
{I'l : a = 1, . . . ,s, L = 1, . . . , oo} such that {j^£} for a fixed a label a partition of Na, in the 
sense that X^L^i '^l^ — ^a- That is, z/£ is the number of times L figures in the partition of 
Na, in other words, the permutation of Na elements is represented by cycles, and vl is the 
number of cycles of length L. The expression for the partition function then is 

Z......=j:F^flf[^(^f , (3.1) 

where the sum is over all conjugation classes V of the group, the double product is the 
contribution from V to the partition function of (A'^i + • • • + A'^s) multispecies bosons, and F-p 
comes from the interaction. By a reasoning completely analogous to that in Refs. 8. 10 . the 
central idea being that the effect of the anyonic "statistics interaction" is to supply a path 
with an extra phase coming from the mutual winding of particles, one concludes that F-p is 
the generating function of winding number probabilities, given by 

Fv=Y^ i^-p({(3jfc})exp(-i7r^aaj.afcQjfc) , (3.2) 

{Qjk} 3<k 

where j or k numbers all the particles, from 1 to J2a^a', o-j is the species that particle j 
belongs to; and P-pdQjk}) is the probability that for a path inducing a permutation from 
class V, the pairwise winding numbers form the set {Qjk}, the distribution of paths being 
the same as for free bosons at inverse temperature f3. 

The winding number Qj^ is defined here as (/>jfc/vr where (j)jk is the angle by which the 
radius vector Tj — r^ rotates. It is an even integer if particles j and k regain their initial 
positions, an odd integer if they get exchanged, but generally a fractional number if one 
or both particles exchange positions with others of the same species. An individual Qjk is 
therefore a continuous variable in general, which seemingly makes the sum in (|3.2() somewhat 
ill-defined. However, it can be observed that 

(i) if particles ji, ■ ■ ■ ,jm make up a cycle (i.e., they are all of species a and get cyclically 
permuted by the path), then it is only the "intracycle winding number" 

m 
p<q=l 

that figures (multiplied by Uaa) in the phase factor in Eq. 1)3. 2|) : 

(ii) if particles ji, . . . ,jm (species a) make up one cycle and /ci, . . . , /c„ (species h) make 
up another cycle, then it is only the "intercycle winding number" 



m n 



Q{jl-jm){k^...k„) = Zm Qjpk, (3-4) 

p=lq=l 



that figures (multiplied by Oab) in the phase factor. 



Now, Qj^...j,n is always an integer (even/odd for m odd/even) and Q[j-^,„j^)[ki...kn) ^^ always 
an even integer. We can then sum over these numbers, with the corresponding probabilities, 
rather than over individual QjkS. Classes of paths characterized by their winding numbers 
form representations of a so-called colored braid group (paths within the same class can be 
continuously deformed one into another). 

We will, as mentioned above, characterize permutations by their corresponding sets of 

1 1 s s 

cycles, and explicitly write out F-p in terms of cycles, as Fj^_^_^"' j^ t^...Lsi...l s) where 
{Lai, ■ ■ ■ , Lau<^} is the set of lengths of cycles pertaining to species a. (The species sequence, 
shown here for simplicity as 1, . . . , s, may in fact miss some members and may not come in 
order.) E.g., -F213 i^ ^^^ generating function for paths involving three particles of species a 
and three particles of species 6, inducing a permutation which interchanges two of the species 
a particles (cycle of length 2) and cyclically permutes all three particles of species h. 

The cluster coefficients are expressed in terms of the F- functions via Eqs. (|2.4|) and (|3.1|) . 
In the thermodynamic limit A — > cx), which we will be interested in, the single-particle 
partition function is Zi{f3) = Aj}? cx /?~"^, so that Zi(L(3) = Zi(f3)/L. Denoting, for brevity, 
Zi = Zi{j3), one gets the following relations: 

X%^ = IgII + \g\ , \%u = G|? ; (3.5) 

•^ ^3 = gG'ni + -G21 + -G3 , A 621 = 2^111 + 4^21 > ^ ^111 = Gill ) (3-6) 

where 

G^ = Ff , Gt\ = [Ff, - l]Zi ; (3.7) 

G^ = F^ Of, = [Ff, - F-^jZi , Gf.'i = [Ff.'i - Ff, - F^ - F^l + 2\Z\ (3.8) 

are the "connected parts" of the -F's, in a sense to be clarified below. 
From the single-species case, it is known [HI E] that 

L-\ 



G2 = n = n(i-T"-) ' (3- 



k=\ 



k 



so, in particular, 



G^ = 1 - 2aaa , G^ = 1 - laaa + ^a^ . (3.10) 

Also, from the exact two-anyon partition function it can be derived that Gii = aaa{(Xaa — !)• 
Since this comes from paths of a single pair of particles without permutation, the same result 
holds for distinguishable particles as well: 

Gt\ = aabiaab - 1) . (3.11) 

Hence, in the second order both the single-species |13| and mixed cluster and virial coefficients 
are known exactly: 

A^62(aii) = -a2(aii) = ^{l - a^f - - , (3.12) 

A^6ii(ai2) = -011(012) = ai2(ai2 - 1) , (3.13) 

Note that 611(012) = 62(012) + 62(1 -1-012), reflecting the fact that the basis wave functions of 
two distinguishable anyons can be chosen to consist of symmetric and antisymmetric functions 
multiplied by the same 012 dependent anyonic prefactor as for identical anyons. 



4 The two-particle approximation 

The physical meaning of the G-functions becomes apparent as soon as they are cast exphcitly 
in terms of winding number probabihties. Consider, e.g., G21, which depends on two statistics 
parameters, an and ai2- (Once we start exphcitly writing down the species-dependent statis- 
tics parameters, we drop the species-related superscripts on the G's as long as this causes no 
ambiguity.) According to the definitions (|3.8|) and H3.2j) ~ H3.4|) . 

G2i(aii,ai2)= Y. ^Qi2Q(i2)(3) (e'"'°"'^^'+"^''5(i^)(«'^ - e'™"'?^^) . (4.1) 

Ql2Q(12)(3) 

Terms with Q(i2)(3) = do not contribute to the sum. In other words, only those paths do 
contribute where the 2-cycle formed by the first two particles, of species 1, actually winds 
around the 1-cycle which is the third particle, of species 2. Generally, a path which is 
disconnected in the sense that there exist two subsets of cycles such that the path involves no 
winding of any cycle in the first subset around any cycle in the second one, will not contribute 
to its corresponding G-function. It is in this sense that we call the latter the connected part 
of the F-function. 

The probability for a path involving two distinguishable particles to have a nonzero wind- 
ing number is inversely proportional to Zi — i.e., it tends to zero in the thermodynamic limit. 
Indeed, a single-particle path, on average, covers an area proportional to A'^, and X'^/A = Z^ 
is the probability of two such paths coming close to each other. This makes it sensible to 
suggest an approximation whose essence is neglecting more than two-particle correlations. 
Specifically, this is an approximation of uncorrelated minimal pairwise windings, in which: 
(i) any intercycle winding is replaced with all possible combinations of pairwise windings be- 
tween one particle belonging to one of the cycles and another one belonging to the other (i.e., 
the fractional pairwise winding numbers are replaced with integers); (ii) only minimally con- 
nected paths are taken into account (e.g., for three particles, if Q12 ^ and Qis 7^ 0, which 
is enough for connectedness, Q23 is assumed to be equal to zero); (iii) windings of different 
pairs are supposed to be mutually uncorrelated [which is, of course, an exact statement with 
respect to pairs (12) and (34), but not with respect to (12) and (13)]. 

Since the generating function for combined probabilities of uncorrelated events factorizes, 
any G-function in this approximation (to be denoted by G) is a sum of products of G^^'s and 
Gil's with different arguments, i.e., a polynomial in the a^b's. In particular, 

G2i(aii,ai2) = 2G2(aii)Gii(ai2) = 2(1 - 2011)012(012 - 1) , (4.2) 

representing a distribution of the 2-cycle windings, G2, and an uncorrelated distribution of 
windings of the third particle around either the first or the second one, 2Gii . 

Contributing to the approximation for Gm are paths with two, but not three of the 
winding numbers nonvanishing; therefore, 

6111(012,013,023) = Gii(oi2)Gii(ai3) +Gii(0i2)Gii(023) + Gii(ai3)Gii(Q23) 

= ai2(ai2 - I)ai3(ai3 - 1) + ai2(ai2 - 1)023(023 - 1) 

+ 013(013 - 1)023(023 - 1) • (4.3) 

One might be tempted to surmise that this expression is in fact exact in the thermodynamic 
limit, the contribution of paths with the third winding number nonvanishing containing an- 
other factor of Z^ . That is of course not the case. If both Q12 and Qis are nonzero. 



particles 2 and 3 are already close to each other, since they are both close to particle 1. Still, 
it is a sensible starting point, as indeed confirmed by the simulations. (For a single species, 
a straightforward diagrammatic interpretation has been adduced JDj. Its generalization to 
many species is possible, as is a certain improvement of the approximation itself |14j.) 

Remarkably, in this polynomial approximation, starting with the third order, all the virial 
coefficients turn out to be statistics independent (in particular, all the mixed virial coefficients 
vanish). Introducing the //-functions which are the deviations from the approximation in 
question, 

Gv = Gv + Hj, , (4.4) 

and using the explicit expressions for the G"s, one gets, in the third order, 

«3(aii) = tjt: - oHiu{aii,au,au) , (4.5) 

021(011,012) = ---^21(011,012) - Hui{au,ai2,ai2) , (4.6) 

0111(012,013,023) = -2i/iii(ai2, 013,023) (4.7) 

(in the first formula, it has been taken into account that i/2i(oii, on) = 0, since for a single 
species, the polynomial approximation for G21 is known E] to be exact), and likewise 
in the higher orders. This is the simplest possible way in which the equation of state can 
interpolate between the Bose and Fermi limits. Absence of corrections of the third and higher 
orders in density is of course a consequence of the multiparticle effects being neglected. The 
corresponding equation of state is characteristic of multispecies exclusion statistics j2]. Thus, 
as well as for a single species [H] , the thermodynamics of exclusion statistics is reproduced by 
anyons as long as only two-particle correlations are taken into account. 

5 Monte Carlo simulation and results 

We now turn to a direct numerical computation of the generating functions, and thence of the 
cluster and virial coefficients. To that end, we Monte Carlo simulate random paths inducing 
a given permutation, with the thermal distribution valid for free bosons. For each path, we 
calculate the intracycle and intercycle winding numbers as defined by Eqs. H3.3() . 1)3. 4|) . The 
resulting winding number probabilities are substituted into the definitions of the generating 
functions. 

We are interested in the thermodynamic limit, in which Z^ — > 0; but in practice, the 
calculations have to be performed at a finite value of Z^ , determined by the choice of /3. 
(Both m and A can be scaled to unity, then Z^ = 27r/3.) A compromise is necessarily 
involved: This value cannot be too big lest finite-size effects emerge, but it cannot be too 
small because, as explained above, the number of useful events tends to zero together with 
it. We have chosen such values (different for different partitions) that typically several paths 
out of a thousand contribute to the connected parts. 

For the (21) partition, 32.1 x 10*5 paths were generated at /? = 0.01 (Zf^ = 0.0628). 
Figure^shows the correction i/2i(oii, 012) to the polynomial approximation G21, Eq. 1)4.2(1 . 
It vanishes, as it should, for 012 = and 012 = 1 (no statistical interaction between the 
two species), as well as for on = 012 (see above). Also, like G21 itself, it is antisymmetric 
with respect to (011,012) 1-^ (1 — oii,l — 012), since the total winding number is always 
odd. The maximum is /72i(0, 0.5) ~ 0.071, which is to be compared with the maximum of 
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Figure 1: i?2i(aii,ai2)- 
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(c) 
Figure 2: The mixed third virial coefficient aiii(ai2, ais, a23) 

"23 = 0, (b) 023 = 0.1, (c) 023 = 0.5. 



-2^111(012,013,023) at (a) 




Figure 3: The mixed third virial coefficient a2i{au,ai2). 



G21 (0,0.5) = 0.500; thus, the relative error of the two-particle approximation is about 15%. 
The absolute error of the MC data themselves, which can be estimated by looking at the 
imaginary part of the simulated function, never surpasses 0.001. 

For the (111) partition, 340 x 10^ paths were generated at /? = 0.03 (Zf ^ = 0.188). We 
illustrate our results with plots of aiii(ai2, ai3,a23) = — 27^111 (ai2,ai3, 023) at fixed values 
of ai3 = 0, 0.1, and 0.5, on Fig.|2l(a), (b), (c), resp. The first of the three cases, where winding 
numbers of one of the pairs do not contribute, reflects the three-particle correlation in the 
windings of the pairs (12) and (13). The extreme value, ffm (0.5, 0.5,0) ~ 0.00138, is to be 
compared with the conditional maximum Gm (0.5, 0.5,0) = 0.06250. The approximation of 
uncorrelated windings per se is, therefore, accurate up to 2%. When all three parameters 
are nonzero, it is not only the above-mentioned correlation but also any paths with all three 
winding numbers nonvanishing that start to contribute to the correction, making the latter 
bigger by almost an order of magnitude. The extreme value is -ffiii(0.5,0.5,0.5) ~ —0.0252 
vs. Gill (0.5, 0.5, 0.5) = 0.1875; the error is 0.00025 at most. Correspondingly, the extreme 
value of the three-species mixed virial coefficient is am (0.5, 0.5, 0.5) ~ 0.0504. 

Our next result. Fig. |31 is the two-species mixed third virial coefficient 021(011,012), 
determined from Eq. (|4.6j) . It includes an antisymmetric [with respect to (011,012) '~^ (1 ~ 
Oil, 1 — 012)] part, H21, and a symmetric one, Hm. The latter being mostly negative, the 
virial coefficient is predominantly positive; however, its extreme negative value, 021(0, 0.5) ~ 
—0.0368 is slightly bigger by magnitude than its extreme positive value, 021(1, 0.5) ~ 0.0343. 

As a final check, we have computed the single-species third-order virial coefficient 03(011), 
and found the result to be consistent with the known formula j2j, 

03(0) = T^ + -^ sin^ vro - (1.65 ± 0.01) x 10"^ sin^ vro , (5.1) 

36 127r^ 



within the error margin of 0.0002. (The magnitude of the error is bigger than that of the sin 
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term, which cannot therefore be seen here.) 

The next logical step is to try to infer the approximate analytic formulas for the mixed 
virial coefficients, either as fits to the MC data or by somehow taking into account multipar- 
ticle correlations in the winding number distributions. It is possible to make some progress 
on both avenues; results will be reported elsewhere |14j . 
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